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SUMMARY 


Anisotropic  plate  theory  was  used  to  calculate  the  anisotropic  stress  concentra¬ 
tion  factors  for  a  composite  plate  (AS/3501-5  graphite/epoxy  composite,  single  ply 
or  laminated)  containing  a  circular  hole.  This  composite  material  was  used  on  the 
highly  maneuverable  aircraft  technology  (HiMAT)  vehicle.  It  was  found  that  the  ani¬ 
sotropic  stress  concentration  factor  could  be  greater  or  less,  than  3  (the  stress  con¬ 
centration  factor  for  isotropic  materials),  and  that  the  locations  of  the  maxinnim 
tangential  stress  points  could  shift  with  the  change  of  fiber  orientation  with 
respect  to  the  loading  axis.  The  effect  of  hole  sire  on  the  stress  concentration 
factor  was  examined  using  the  Point  Stress  criterion  and  the  Averaged  Stress  Cri¬ 
terion.  The  predicted  stress  concentration  factors  based  on  the  two  theories  com¬ 
pared  fairly  well  with  the  measured  values  for  the  hole  size  0.3175  cm  (1/8  in).  It 
was  also  found  that  through  the  lamination  process,  the  stress  concentration  factor 
could  be  reduced  drastically,  indicating  an  improvement  in  structural  performance. 


INTRODUCTION 

When  an  infinite  isotropic  plate  containing  a  circular  hole  is  subjected  to 
remote  uniaxial  tension,  the  tangential  stress  along  the  boundary  of  the  circular 
hole  will  reach  a  value  three  times  the  remote  tensile  stress  at  two  boundary  points 
lying  on  the  hole  diameter  perpendicular  to  the  loading  axis.  Namely,  the  maximum 
tangential  stress  concentration  factor  is  3,  %diich  is  independent  of  the  hole  size. 
For  anisotropic  materials,  such  as  fiber-reinforced  composite  materials,  the  picture 
is  entirely  different.  The  value  of  the  maximum  tangential  stress  concentration  fac¬ 
tor  for  a  composite  plate  can  be  greater  or  less  than  3,  and  the  locations  of  the 
maximum  stress  points  could  shift  depending  on  the  loading  direction  and  the  fiber 
orientations • 

This  report  calculates  the  tangential  stress  distribution  around  a  circular  hole 
in  a  composite  plate  (single  ply  or  laiainated)  and  examines  how  the  maximum  stress 
points  shift  with  the  fiber  orientations.  In  addtion,  the  dependency  of  the  stress 
concentration  factor  on  the  hole  size  is  examined.  In  the  analysis,  the  coa^site 
system  is  treated  as  a  continuous  anisotropic  plate  having  effective  elastic  proper¬ 
ties.  The  calculated  stress  concentration  factors  are  then  compared  with  experimen¬ 
tal  data. 


NOMENCLATURE 

small  fixed  distance  ahead  of  the  hole  boundary  used  in  the  Average. 
Stress  Criterion 

characteristic  distance  ahead  of  the  hole  boundary  used  in  the  Point 
Stress  Criterion 

modulus  of  elasticity  of  single  ply  in  aucis-l  direction  or  modulus  of 
elasticity  of  anisotropic  plate  in  axis-1  direction 
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modulus  of  elasticity  of  single  ply  in  axi8'*2  direction  or  modulus  of 
elasticity  of  anisotropic  plate  in  axi8-2  direction 

modulus  of  elasticity  of  laminated  composite  plate  in  axis-l  direction 

modulus  of  elasticity  of  laminated  composite  plate  in  axi8-2  direction 

modulus  of  elasticity  of  anisotropic  plate  in  a  direction 

modulus  of  elasticity  of  single  ply  in  fiber  direction 

modulus*  of  elasticity  of  single  ply  transverse  to  fiber  direction 

shear  modulus  of  single  ply  associated  with  {l#2}  system 

shear  modulus  of  single  ply  associated  with  {L,T}  system 

stress  concentration  factor 

stress  concentration  factor  at  a  «  v/2 


coordinate  axis  in  fiber  direction 

number  of  composite  layers 

radius  of  circular  hole 

rectangular  Cartesian  coordinates 

coordinate  axis  transverse  to  fiber  direction 

rectangular  coordinate  system 

angular  cooordinate 

angle  between  axis  1  and  axis  L 

complex  roots  of  the  anisotropic  plate  characteristic  equation 

Poisson *8  ratios  of  single  ply  with  respect  to  {lf2}  syst^ 

Poisson's  ratio  of  single  ply  with  respect  to  {L,t}  system 

Poisson's  ratios  of  laminated  composite  plate  with  respect  to  {l#2} 
system 
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The  HlMAT  cwnposlte  plate  (A5/3501-5  graphite/epoxy  coii?>08ite)  is  made  up  o£  34  plies 
(D  m  34,  ply  thickness  *  0.0133  cm  (0.00525  in))  with  the  following  fiber 
orientations: 


14  plies  of  6  ■  +50*  fiber  orientation 
14  plies  of  9  *  -50*  fiber  orientation 
6  plies  of  0  ■  +35*  fiber  orientation 


The  ply  engineering  elastic  constants  with  respect  to  the  principal  elastic  axes 
{L,t}  are  given  by 


Ei,  -  137.90  GPn  (20  x  10*  Ib/in^) 
E^  -  10.27  GPa  (1.49  x  10*  Ib/in^) 
Glt  «  2.41  GPa  (0.35  x  10*  Ib/in^) 
“  0.3 


(11) 

(12) 

(13) 

(14) 


Using  equations  (11)  to  (14),  the  ply  elastic  constants  with  respect  to  axes  {1,2} 
can  be  calculated  from  equations  (1)  to  (5).  Por  0  ■  ±50*  filser  orientations. 


B^(i50*)  «  7.45  GPa  (1.0809  x  10*  Ib/in^) 

E2(i:50*)  ■  8.43  GPa  (1.2237  x  10*  Ib/in^) 

Gi2(*50)  -  8.46  GPa  (1.2274  x  10*  Ib/in^) 

v,2(±50)  -  0.5676 

For  0  •  +35*  fiber  orientation, 

Ei(+35*)  -  9.59  GPa  (1.3904  x  10*  Ib/in^) 

B2(+35*)  a  7.40  GPa  (1.0732  x  10*  Ib/in^) 

6|2(+35*)  ->  6.91  GPa  (1.0025  x  10*  Ib/in^) 


(15) 

(16) 

(17) 

(18) 


(19) 

(20) 
(21) 


*Hote  that  the  mixture  rule  does  not  give  the  relationship  v^2E2  “  ''21*1  * 


Vt2(+35»)  -  0.6671 


(22) 


Using  equations  (15)  to  (18)  and  equations  (19)  to  (22),  the  engineering  elastic 

constants  for  the  HiMAT  coii5)osite  plate  with  respect  to  {1,2}  system,  can  be  calcu¬ 
lated  from  equations  (6)  to  (9).  j  >  ^ 


*  n  =1  <+50*)  +  ^  E,  (-50«)  +  E,  (+35*) 

■  7.83  GPa  (1.1356  x  io®  Ib/in^) 

®2  “  ^  E2  (+50»)  +  ^  E2  (-50«)  +  £3  (+35* ) 

•  8.25  GPa  (1.1972  x  10^  Ib/ln^) 

*^12  -  ^  Gi2  (+50*)  +  ^  Gi2  (-50»)  +  G,2  (+35*) 

“  8.19  GPa  (1.1877  X  io6  ib/in2) 

^12  -  n  '’12  (+50«)  +  ^  v,2  (-50«)  +  v^2  (+35»)  -  0.5851 


(23) 


(24) 


(25) 

(26) 


TANGENTIAL  STRESS  AROUND  A  CIRCULAR  HOLE 


For  an  anisotropic  plate  containing  a  circular  hole  subjected  to  remote  uniaxial 
tensile  ctress  a^,  acting  at  an  angle  ^  with  respect  to  the  principal  elastic  axis  1 

of  the  plate  (fig.  2),  the  tangential  stress,  Oo  (or  tangential  stress  concentration 
tactor,  K  =  Ofj/aJ  along  the  circular  hole  boundary  may  be  expressed  as  (ref.  3, 
p.  174) 


^  ^  ♦  +  (l«+n)  sin2  k  cos2  a 

00  1 

+  [(1+n)  cos2  ^  -  k  8in2  sin2  a 

+  k  +  n)  sin  cos  iP  sin  (X  cos  q| 

where  E^  is  the  modulus  of  elasticity  in  the  o  direction  (fig.  2)  given  by 


(27) 


(28) 


where  k  anc2  n  are  defined  by 


n 


(30) 


»  -1(1.,  ♦  ^ 

where  i  =  and  and  U2  are  the  complex  roots  of  the  anisotropic  plate  charac¬ 

teristic  equation 

+  ij  »  0  <31) 

For  isotropic  materials  k  *  1  and  n  «  2,  and  the  stress  concentration  factor  K 
(eq.  (27))  reduces  to 


''a 

K  “  —  *  1-2  cos  2(a  -  ^)  (32) 

which  gives  K  *  -1  at  (a  -  ♦)  «  0  or  w,  and  K  »  3  at  ( o  -  ^)  «  ±  ir/2 


HOLE  SIZE  EFFECT 


Consider  a  circular  hole  of  radius  R  in  an  infinite  anisotropic  plate  as  shown 
in  figure  3-  If  the  remote  uniform  stress  is  applied  in  the  y-axis  (or  axis  1) 

direction,  then  the  normal  stress  Oy  in  the  y-axis  direction  at  a  point  on  the  x-axis 
(or  axis  2)  in  front  of  the  hole  may  be  approximated  by  (ref.  4) 

^,..0,  .  I--  j,  .  (B).  .  ,{£)* 

-((1  +  n)  - 
(x  >  R) 


At  the  hole  boundary  x  »  R  ,  equation  (33)  gives  the  stress  concentrations  factor 

V2 


a„(R,0) 

KV2  ' 


0F=Tr/2,  <fr=0 


1  +  n 


(34) 


whicr  agrees  with  equation  (27)  by  setting  a  =  Tf/2  and  -  0.  Equation  (34)  gives  a 
constant  value  of  K7r/2  the  same  material  regardless  of  the  hole  size.  As  the 

hole  size  decreases,  the  stress  concentration  factor  K^/2  decreases  and  finally 

^^ity  (that  IS,  a  plate  without  a  hole).  For  composite  materials,  the 
hole-size  effect  on  the  stress  concentration  factor  Kn/2  becomes  significant  when 
the  hole  diameter  becomes  less  than  3.048  cm  (1.2  in)  (ref.  4). 

To  account  for  the  hole-size  effect  in  equation  (34),  two  failure  criteria  were 
advanced  by  Whitney  and  Hulsmer  (ref.  5)  for  coi'q>oslte  materials.  Two  theories  are 
described  in  the  following  sections. 


6 


Point  Stress  Criterion 

The  point  stress  criterion  assumes  that  the  failure  will  occur  when  the  stress 
Oy  {x,0)  at  a  certain  small  fixed  distance  (or  characteristic  distance)  dQ  ahead  of 

the  hole  boundary  first  reaches  the  tensile  strength  Of  of  the  material  (or  tensile 
strength  of  the  plate  without  a  hole)  (Fig.  3)*  Namely, 

I 

Using  this  criterion  and  equation  (33),  the  stress  concentration  factor, 
can  be  written  as  (ref*  5)  • 

^  2  34  (1+n)-3/  ^ 

where 

=5-^5;  <37) 

For  a  large  hole  (that  is  Ci+D,  equation  (36)  gives 

I  -  Kt/2  -  1  +  n  (38) 

As  the  hole  size  decreases  (that  is,  Ci+O),  equation  (36)  reduces  to 


which  corresponds  to  the  case  of  no  holee 


Average  Stress  Criterion 


The  average  stress  criterion  assumes  that  failure  will  occur  when  the  average 
value  of  (^(x,0)  over  some  small  fixed  distance  a0  ahead  of  the  hole  boundary  first 
reaches  the  tensile  strength  of  the  material  (without  a  hole)  (fig.  3).  Namely, 


(x,0)  dx 


(40) 


Using  equations  (33)  and  (40),  the  stress  concentration  factor  Kwo  can  be  written  as 
(ref.  5) 
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(41) 


(2) 

Kit/2 


2(1  -  52) 
(1  +  52) 


^2  -  52  -  4  +  +  n)  -  3]  (5^  -  52)]* 

^2  +  5^+  t(1  +  n)  -  3)52} 


(42) 


where 

^2  ■  R  +  ao 

For  a  large  hole  (that  is,  52^1 )*  equation  (42)  reduces  to  equation  (34). 

(2)  I 

K^/2  j  ■  Kir/2  ”  1  +  n 

And  for  the  case  without  a  hole  (that  is#  52^^)^  equation  (42)  gives 


(43) 
Namely, 

(44) 


(45) 


The  values  of  the  characteristic  distances  d^  and  a^  are  determined  by  means  of  curve 
fitting  of  the  experimental  data  obtained  from  tensile  tests  of  rectangular  specimens 
containing  holes  of  different  sizes.  Fr<»  literature  surveys  (refs.  4  and  5),  the 

values  of  do  and  Sq  for  con^site  materials  are  in  the  ranges  of 


do  •  0.076  to  0.127  cm  (0.03  to  0.05  in) 
ao  **  0.381  cm  (0.15  in) 


♦ 


EXFERIMEHTS 


The  effect  of  hole  size  on  the  stress  concentration  factors  predicted  from 
the  two  theories  discussed  previously  were  examined  by  performing  simple  coupon 
tests.  The  widths  W  of  the  two  rectangular  specimens  were  W  »  3.81  cm  (1.5  in) 
and  W  «  1.905  cm  (0.75  in).  Both  of  the  specimens  contain  small  circular  holes  of 
2R  -  0.3175  cm  (1/8  in).  The  specimens  were  tested  in  a  tensile  test  machine  up  to 
failure  under  uniaxial  tension  in  the  axrs-1  direction,  and  the  remote  stress  at 

the  time  of  failure  was  calculated  for  each  specimen.  One  specimen  without  a  hole 
was  also  loaded  up  to  failure  in  a  similar  manner  to  obtain  the  tensile  strength.  Of, 
of  the  material.  The  stress  concentration  factor  Ktr/2  each  specimen  with  a  hole 
was  calculated  from 


RESULTS 


By  replacing  Ej,  ''12I'  respectively,  with  {i^,  E2t  G\2>  '’12! 

equations  (27)  and  (28),  the  tangential  stresses  around  a  circular  hole  in  a 
laminated  HiMAT  composite  plate  were  calculated  for  three  loading  cases:  ^  »  0 
(loading  in  axis-1  direction),  ♦  -  »/4,  and  ♦  =■  «/2  (loading  in  axis-2  direction). 
Figure  4  shows  the  plots  of  a  a  for  both  the  laminated  HiMAT  composite  (eq.  (27))  and 
isotropic  materials  (eq.  (32))  when  the  plate  is  under  uniaxial  tension  in  the  com¬ 
posite  elastic  axis-l  direction.  The  maximum  stress  concentration  factor  K  for  the 
laminated  HiMAT  coii?>osite  reached  the  value  2.43  (less  than  3)  at  f « ur  locations 
(a  »  i  65®  and  o  «•  r  115®)  instead  of  two  locations  (o  *  i  */2)  for  he  isotropic 
case.  When  the  loading  axis  is  7  -  ir/4  oblique  to  the  composite  axis  1  (fig.  5),  the 
stress  concentration  factor  K  reaches  the  value  of  3.49  (greater  than  3)  at  two  loca¬ 
tions  (a  =  135®  and  o  »  -45®).  When  the  loading  axis  is  parallel  to  the  conqposite 
elastic  axis  2  (fig.  6),  the  stress  concentration  factor  K  reaches  the  peak  value  of 
2.44  (less  than  3)  at  four  locations  (a*  *25®  and  a  =  *155®). 

For  comparison  purposes,  similar  calculations  were  made  for  a  single  ply  of  the 
HiMAT  conmosite  using  equation  (27).  When  loading  is  along  the  fiber  direction 
(axis  1,  fig.  7),  the  stress  concentration  factor  K  reaches  10.20  at  four  locations 
(a  -  *80®  and  a  »  *100®).  When  the  loading  is  45®  oblique  to  the  fiber  direction 
(fig.  8),  K  reaches  as  high  as  27.78  at  two  points  (a  -  114®  and  a  -  -66®).  When  the 
loading  direction  is  transverse  to  the  fiber  direction  (fig.  9),  the  value  of  K 
reaches  6.81  at  four  points  (u  ■  -55®  and  o  ■  *125®).  Tables  1  and  2  summarise  the 
above  results  for  both  laminated  and  single-ply  HiMAT  composites.  Notice  that  the 
laminated  coii?)08ite  has  much  lower  stress  concentration  factors  than  the  single  ply, 
indicating  the  in^rovement  of  the  structural  performance  through  the  lamination 
process. 

Figure  10  shows  the  plots  of  equations  (36)  and  (42)  for  the  laminated  HiMAT  com- 


posite 


n  »  j 

N3  '  ~ 

1 

Cl 

V=2  / 

■)-ti 


taking  the  characteristic  distances  6^  *  0.1016  cm 


ini  F  —  V.  - — ^ 

a  hole  diameter  of  0.3175  cm  (1/8  in),  the  predicted  values  of  Kit/2  are  summarized 
below. 


cm  (in) 

0.0106  (0.04) 
0.1270  (0.05) 


»o* 

cm  (in) 


0.3810  (0.15) 


1.4367 

1.3152 

1.3498 


The  stress  concentration  factors  Kr/2  determined  from  the  uniaxial  tensile  tests 
(loading  in  axis-1  direction)  of  the  laminated  HiMAT  conq>osite  specimens  are  tabu- 
lated  below. 
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Speclnan 


2R, 

cm  (in) 


1 

2 


W, 

cm  (in) 


KV2 


3.810  (0.75) 
1.905  (1.5) 


0.3175  (1/8) 
0.3175  (1/8) 


1.333 

1.220 


The  above  two  data  points  are  plotted  in  figure  10  and  compare  fairly  well  with 
the  predicted  curves,  especially  with  the  curves  d^  -  0,127  cm  (0.05  in)  and  a^  = 
0.381  cm  (0.15  in). 


CONCULDING  REMARKS 


Anisotropic  plate  theory  was  used  to  calculate  the  anisotropic  stress  con¬ 
centration  factors  for  single-ply  and  laminated  HiMAT  (hi^ly  maneuverable  aircraft 
technology)  research  aircraft  composite  plates,  each  of  which  contained  a  circular 
hole.  The  analysis  showed  that  the  anisotropic  stress  concentration  factor  could  be 
greater  or  less  than  three  for  anisotropic  materials,  (three  is  the  stress  concentra¬ 
tion  factor  for  isotropic  materials)  and  the  locations  of  the  maximum  tangential 
stress  points  could  shift  with  the  chan^  of  fiber  orientation  in  relation  to  the 
loading  axis.  It  was  found  that  throu^  the  lamination  process  the  stress  concentra¬ 
tion  factor  could  be  reAiced  drastically,  and  thereby  the  structural  performance 
could  be  improved.  The  effect  of  the  hole  size  on  the  stress  concentration  factor 
was  examined  in  the  light  of  the  Point  Stress  Criterion  and  the  Averaged  Stress  Cri¬ 
terion.  The  stress  concentration  factors  calculated  from  these  theories  agree  fairly 
well  with  the  measured  values  for  the  hole  size  0.3175  cm  (1/8  in). 
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TABLE  1 .  -  STRESS  CONCENTRATION  FACTORS  FOR  A  LAMINATED  HiMAT  COMPOSITE  PLATE, 
34  PLIES  (t±50*Ji4,  [-50*J^4,  t+35*]g) 

Orientation  of  loading  axis,  ^ .  0*  45*  90* 

Stress  concentration  factor,  2*43  3«49  2.44 

Locations  of  peak  tensile  stress,  a.*  ±65*,  ±115*  +135*,  -45*  ±25*,  ±155* 

TABLE  2.  -  STRESS  CONCENTRATION  FACTORS  FOR  A  SINGLE  PLY  OP  HiMAT  COMPOSITE 

Orientation  of  loading  axis,  .  0*  45*  90* 

Stress  concentration  factor,  K  10.20  27.78  6.81  i 


Locations  of  peak  tensile  stress,  a,. I  ±80*,  ±100*  +114*,  -66*  I  ±55*,  ±125’ 


Figure  !•  Single  ooa-  Figure  2m  Tension  et  an  angle  to  a 
posite  ply  under  oc^ined  principal  elastic  axis  1  of  an  ani- 
loadingm  sotropic  plate  9fitb  a  circular  hole 


Figure  Tangential  stress 


Figure  7 .  Tangential  stress  distri^ 
but  ion  around  a  circular  bole  in  a 
single^ply  composite  plate*  ♦  •  • 


distraction  around  a  circular 
hole  in  a  single~ply  composite 
plate*  ♦  *  45* * 


Figure  9.  Taageatial  stress 
dlstribatioa  aroaed  a  circular 
hole  in  a  single-pig  coaposite 
plate,  ♦  •  9(F, 
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